Diffraction-free nonevanescent nano-beams using the Presnel-waveguide concept 
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Strong intensity attenuation limits the use of conventional diffraction-free optical elements. We 
show a possible solution to the exponential intensity attenuation limiting the use of Fresnel-type 
diffraction-free nanometer-scale optics by using materials with appropriately chosen refractive index. 
Such beams may be applied for technical and physical problems. 
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I. INTRODUCTION 

Resolving micro- and nanoscopic objects is a challeng- 
ing problem and it has been an important field of research 
in the last decades. Although a simple microscope is 
sufficient to magnify microscopic images, its resolution 
capability is confined by the Abbe diffraction limit: 



d = 



A 



2n sin(0) : 



(1) 



where A is the wavelength, n is the refractive index while 
the angle of incidence is denoted by 9. Therefore diffrac- 
tion is an ultimate barrier in traditional microscopy. 

Many solutions have been proposed to overcome this 
limit, that is, to eliminate its effect, such as photoactivat- 
able localization microscopy [l| (PALM), stimulated emis- 
sion depletion]!] (STED), surface plasmon polaritonsQ 
(SPPs), various nanostructures such as the recently dis- 
covered microspheres or the diffraction- free wave solu- 
tions p-[lJL Ultrashort diffraction- free pulses can also be 
produced jl8l. l20l l2ll | . Diffraction- free light beams can be 
defined as beams with no diffractive broadening through- 
out the free-space propagation. This can be expressed as 
follows: for all a;, y and z coordinates the average trans- 
verse components It — St of the Poynting vector St of 
the beam must fulfill the following condition fl5| : 



V-I-i 







(2) 



In [18| we described a diffraction-free arrangement that 
consists of an array of slits illuminated by TM- waves with 
periodically altering phase; however, it was also pointed 
out that this diffraction-free feature implies an exponen- 
tial decrease in the intensity, i.e. behind the slit array 
evanescent beams can be observed. 

Evanescent beams are described through an exponen- 
tial term: 



U(x,y, z) = F(x,y, z) ■ exp(ifiz), 



(3) 



where z is the coordinate along the propagation direc- 
tion, and (3 is the wavevector component in the direction 
of propagation. If j3 is purely imaginary, the beam be- 
comes evanescent; however, if this wavevector component 
remains real, propagating waves are generated. 



This was the basic idea of Ruschin and Leizer |16T ]. 
In their paper they concentrated on Bessel-type beams 
which arc generally treated as propagating waves. They 
focused on the problem of generating evanescent Bcssel- 
beams, therefore in their work they used electromagnetic 
solutions of the form 

U m (r, 0, z) = exp(i^m)J m (ar) exp(-/3'z), (4) 



where /3' = a 2 — n 2 k 2 is the propagation constant, 
(r, (j), z) are cylindrical coordinates, J m is the Bcsscl func- 
tion. Whenever j3' is real, i.e. a 2 > n 2 k 2 , this func- 
tion decreases exponentially; however, if /3' is imaginary, 
the exponential part describes a propagating component. 
This can be achieved for example by changing the refrac- 
tive index of the materials. 

Another important result in this field was published by 
Yacob Ben-Aryeh[T3l- This paper analyzes theoretically 
the microsphcre-based field conversion from evanescent 
to propagating waves Q ■ Information on the nanostruc- 
ture of a nano-corrugated metallic thin film is included 
in the evanescent waves produced by the plane EM waves 
transmitted through the structure. This information can 
be extracted from the evanescent wave by microspheres 
located above the surface that collect and convert them 
into propagating waves. The conversion effect is also the 
consequence of the refractive index difference of the ma- 
terials. 

Nondiffracting beams have some special features which 
may be used in technical as well as physical applications 
such as imaging, particle acceleration or light motors [231- 
[28| . In these applications, however, the use of nonevanes- 
cent beams would be preferred. 

In our paper we show that the production of prop- 
agating (nonevanescent) diffraction-free nano-beams us- 
ing the Fresnel- waveguide concept is feasible. Our geom- 
etry is similar to a simple grating; however, the phase 
difference between the EM fields falling to the adjacent 
sub-wavelength slits and the near-field treatment of the 
problem distinguish our case from that of a simple grat- 
ing. The paper is arranged as follows. In Section |TT] we 
present the basic idea and the results of our theoretical 
analysis. Section IlIII contains the results of our numeri- 
cal simulations using a finite element analysis code. We 
present different geometries and conclude that critical re- 
fractive indices separating the evanescent and propagat- 
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ing regions depend on the transverse characteristic sizes 
of the geometries. Finally, in Section HVl conclusions are 
drawn and some final remarks are made. 



II. THEORY 

In [l8| we showed that it is possible to produce 
diffraction-free nano-beams in the subwavelength regime 
using the Fresnel- waveguide concept. Throughout our 
analysis we modeled an array of one-dimensional slits 
(Fig. [IJ of width 2a in a screen of thickness b illuminated 
by TM-polarized plane waves of wavelength A. The phase 
4> m of the wave falling on the m-ih slit of the array was 
adjusted as <j> m — m ■ it, m — 0, ±1, ±2 . . . These phase 
jumps are very important as the phase difference of the 
fields distinguishes our problem from that of a grating. 
In our earlier paper we observed that the width of the 
central beam remains the same for great distances and 
the beam shape is sinusoid in the direction perpendicular 
to the propagation direction. However, we also presented 
that the intensity decreases exponentially. In this section 
we find a solution to this problem. 

The beam leaving the slit array propagates in homo- 
geneous material; in this case the propagation can be de- 
scribed by the Hclmholtz equation and we may use the 
scalar field U (r, t): 



{V 2 + k 2 )U{r 1 t) = 0, 



(5) 



where V 2 is the Laplacian and k is the wavenumber. 
Owing to the geometry of the structure and the TM- 
polarization (z-invariance) of the incident light, in our 
case there is no z-dependency. We also disregard time de- 
pendency. Moreover, U(x,y) can be split into two func- 
tions, F(x) and G(y), i.e. variables can be separated. 

Now we assume an exponential form for the y- 
dependent part: G(y) ~ exp(i/c y y), where k y is the 
y-component of the wavevector (propagation constant). 
Substituting this into ([5]) we get: 



dx 2 



F(x) = -(k 2 - k 2 y )F(x). 



(6) 



Considering that a harmonic function is expected in the 

x direction: Fix) ~ sin(^//c 2 — k 2 x), by solving ((6]) for 

a beam with harmonic shape in the transverse direction 
and exponential function in the longitudinal direction, 
one may obtain the following function: 

/ 27ra;\ 

U(x, y) = F(x) ■ G(y) ~ sin f — j ■ exp(ik v y), (7) 



where k„ 



2 ^-) 2 , and Aq is the wavelength 



in vacuum. From this latter formula one may draw the 
conclusion that there are some situations (specific geo- 
metric dimensions and wavelengths) in which the oth- 
erwise evanescent waves become propagating, namely if 



k 2 < the beam is evanescent while in the case of k 2 > 
a propagating beam can be obtained. This leads to the 
following inequalities (note that all the quantities - n, a, 
Ao - are positive real values): 

• evanescent beams are described by: 



n < 



Ao 
4a 



• propagating beams are described by: 

An 



n > 



la 



We also note that this formula was derived for the ge- 
ometry in Figure Q] analyzed in our earlier paper [l8[. A 
more general formula can be written through the trans- 
verse characteristic dimension of the geometry, e.g. in our 

case A = ( — J , therefore 
V 4a / 




(8) 



This deduction leads to the conclusion that in our ge- 
ometry it is possible to achieve diffraction-free and prop- 
agating waves by 1) increasing the refractive index 2) 
decreasing the wavelength or 3) increasing the transverse 
characteristic dimension of the geometry. In the follow- 
ing section we show some simulation results for various 
refractive indices and two different geometries (transverse 
characteritic sizes). 



III. NUMERICAL SIMULATIONS 

In order to verify our result, we performed finite ele- 
ment analysis simulations for various geometries (trans- 
verse characteristic dimensions) and refractive indices. In 
the followings we present some results of our simulations. 

The first simulated geometry was the same we analyzed 
in Section|TT](Fig.[T]). In this case we used an array of slits 
of width 2a in a screen of thickness b. It is illuminated 
from region I by TM polarized field that is perpendicular 
to the array of slits, and the phase of which was changed 
by 7r in the adjacent slits, i.e. <p m = m ■ ix. In this case 
no distance was assumed between the individual slits, 
therefore this geometry is quite hypothetic. We will show 
a more realistic model hereinafter. 

The aim of this paper was to find a solution to the ex- 
ponential decrease of the wave leaving the slit (in region 
III). To prove that the solution proposed in Section HT1 
leads to nonevanescent (propagating) waves we plotted 
the energy flow S y along the y axis behind the screen at 
x = for various refractive indices (Fig. [2]). For compar- 
ison the same quantity for a real material (germanium, 
n = 4.7096 @ A = 800 nm @) was also plotted. In 
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FIG. 1: The original geometry used in In the present 

article we used the following values: wavelength in vacuum 
Ao = 800 nm, slit width 2a — 80 nm, screen thickness 
b = 50 nm. The "individual" slits are separated by perfectly 
conducting material. The phase of the EM field illuminating 
the slits changes periodically (in the figure it is indicated by 
the plus and minus characters) 



this case the characteristic dimension of the geometry is 

A = | — V therefore the critical refractive index of the 

geometry that separates the evanescent and propagating 
case can be formulated as n = Considering a wave- 
length of A = 800 nm and slit width 2a = 80 nm, the 
critical refractive index separating the evanescent and 
propagating cases is n = 5. Indeed, Figure [5] confirms 
this value. Note that the transition from evanescent to 
propagating waves is accompanied by an oscillating ef- 
fect, i.e. in the case of materials with refractive indices 
close to the critical value the y-component of the Poynt- 
ing vector S y of the wave leaving the slit array oscillates 
while decreases; this is shown for n — 4.9 in the inset of 
Fig.[2J The effect is not so apparent for germanium partly 
because it is further from the critical value and partly be- 
cause the material has nonzero conductivity (c = 217). 

To verify that the beam preserves its diffraction-free 
quality, in Fig. [3] we plotted the energy flow S y along the 
x direction at y = 32a for n = 4.9 (evanescent region) 
and n = 4.99 (nonevanescent region). We only present an 
example of the beam shape at y = 32a, our calculations 
showed, however, that the energy flow S y of the beam 
preserves its shape in the region \y\ > 2a, i.e. this phe- 
nomenon cannot be attributed to self-imaging. Although 
separate y axes must be assigned to the two cases because 
of the magnitude difference, the diffraction-free behavior 
can be observed in both cases; however, in the propagat- 
ing case the originally single peak becomes double which 
shows a better localization. This phenomenon, however, 
cannot be observed in the E x or H z components from 
which S y is calculated. The only difference that can be 
observed - apart from the intensity increase - is a change 
in the sign of the E x field. In the case of low refractive 
indices the sign of E x and H z are opposite while in the 



high-refractive-index case their sign is the same which 
results in a double peak. 

In Figure U we present a more realistic geometry. In 
this case the separation s between the slits is nonzero. 
Actually, the separation can be arbitrary, and still the 
diffraction-free quality remains; as an example, in our 
simulations we used s = 4a. The increase of the sep- 
aration affects the characteristic dimension: it becomes 
A = 4a+2s . The width of the diffraction- free beams also 
changes: w — 2a + s. In fact, the former value is the 
period of the slit system; the latter value is presented 
in Figure [5l In the figure we show the energy flow S y 
along the x direction at y = 32a for n = 1.5 (evanes- 
cent region) and n = 1.66 (nonevanescent region). As 
it is expected from Eq. ((?]), the numerical model showed 
that in order to produce propagating beams, increasing 
the characteristic dimension two times allows us to de- 
crease the refractive index two times. It can be seen that, 
similarly to Fig. [31 separate y axes must be assigned to 
the two cases because of the magnitude difference, the 
diffraction-free quality of the beam still can be observed, 
and better localization can be observed for the propa- 
gating case. However, the width of the beam is much 
greater (w = 2a + s) than in Fig. [3J and the critical re- 
fractive index is reduced to n « 5/3 which corresponds 
to a characteristic dimension of A = ia + 2s . 

For comparison we also show the energy flow S y along 
the y axis at a; = for various refractive indices (Fig. . 
In this case it is not only possible to find a real material 
with appropriate refractive index but these materials arc 
liquid therefore they can be used in real experiments. As 
an example, in the simulation we used carbon disulphidc 
(n = 1.6295, a rts @ A = 800 nm). It is clear from the 
inset of Fig.[6]that while S y converges to zero for n = 1.5, 
the value to which S y converges for CS2 differs from zero, 
i.e. it is possible to realize diffraction-free nonevanescent 
beams with this geometry. 

As a final remark we note that our results are in ac- 
cordance with the Abbe-limit (|T|). If sin (9 = 1, the ex- 
pression changes to d = X/(2n). Substituting the values 
used for the simulations, we find that the Abbe-limit for 
11 = 5/3 is d = 240 nm, and for n = 5 it is d = 80 nm. 
These are exactly the transverse characteristic sizes we 
used for the two cases. Consequently, there is a clear link 
between the evanescence limit and the Abbe-limit. 

Our results showed that varying the refractive index 
of the material in which the diffraction-free beams are 
generated leads to nonevanescent beams. These non- 
diffracting beams show some unique properties useful for 
both physical and technical applications. Nondiffracting 
beams may be found applications in imaging as it wa ver- 
ified that extremely long focal depths can be produced in 
imaging realized with such beams [Hj]. Particle accelera- 
tion is another field of research in which diffraction-free 
beams may be used; these beams can serve as alternative 
accelerating fields owing to the strong longitudinal com- 
ponent of the electric field [26[ . It is also possible to ma- 
nipulate ensembles of particles in multiple planes [13, HU 
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FIG. 2: Energy flow on the propagation axis of the central slit calculated behind the slit array shown in Fig. [T] for various 
refractive indices. Slit width 2a = 80 nm, wavelength in vacuum Ao = 800 nm, screen thickness b — 50 nm. For comparison we 
plotted the same quantity for a real material with high refractive index (germanium, n ~ 4.7096). 



which can be realized due to the self-reconstruction abil- 
ity of the nondiffracting beams. As a final example, we 
draw attention to the possibility of realizing light motors 
driven by the diffraction-free beam. However, extracting 
this nondiffracting nonevanescent beam from the higher 
refractive index material to air or any lower refractive 
index one is a challenging problem. 

IV. CONCLUSION 

In our paper we show that the production of prop- 
agating (nonevanescent) diffraction-free nano-beams us- 
ing the Fresnel-waveguide concept is feasible. The an- 
alytical expression derived for nonevanescent beams in 
Section [TlJ similarly to the well-known Abbe-limit, con- 



nects the wavelength in vacuum, the refractive index and 
the characteristic size of the slit system. We presented 
the results of finite element simulations for two slightly 
different geometries and various refractive indices as well 
as real materials. Our numerical results support the find- 
ings presented in Section [Til 
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FIG. 3: Energy flow S y along the x direction at y — 32a for refractive indices n — 4.9 and n — 4.99. Parameters: wavelength 
Ao = 800 nm, slit width 2a — Ao/10, screen thickness b = 50 nm. Transverse characteristic dimension of the geometry is 
A = Note that, although this is only an example of the beam shape at y = 32a, our calculations showed, however, 

that the energy flow S y of the beam preserves its shape in the region \y\ > 2a in both cases, i.e. this phenomenon cannot be 
attributed to self-imaging. 
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FIG. 5: Energy flow S y along the x direction at y — 32a for n — 1.5 (evanescent region) and n = 1.66 (nonevanescent region). 
Parameters: wavelength Ao = 800 nm, slit width 2a — Ao/10, screen thickness b = 50 nm, separation s = 4a. Characteristic 
dimension: A = 2 + , critical refractive index n ss 5/3. 
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FIG. 6: Energy flow S y along the y axis at x = for refractive indices n = 1.5, 1.64, 1.65. Parameters: wavelength Ao = 800 nm, 
slit width 2a = Ao/10, screen thickness b = 50 nm, separation s = 4a. Characteristic dimension: A = 2 "t , critical refractive 
index n « 5/3. For comparison we plotted the same quantity for a real material with moderately high refractive index (carbon 
disulphide, n ~ 1.6295). 



